Sodium channel selectivity and conduction: Prokaryotes have devised 
their own molecular strategy 
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Striking structural differences between voltage-gated sodium (Nav) channels from prokaryotes (homotetramers) 
and eukaryotes (asymmetric, four-domain proteins) suggest the likelihood of different molecular mechanisms for 
common functions. For these two channel families, our data show similar selectivity sequences among alkali cations 
(relative permeability, Pion/PNa) and asymmetric, bi-ionic reversal potentials when the Na/K gradient is reversed. 
We performed coordinated experimental and computational studies, respectively, on the prokaryotic Nav chan- 
nels NaChBac and NavAb. NaChBac shows an "anomalous," nonmonotonic mole-fraction dependence in the pres- 
ence of certain sodium-potassium mixtures; to our knowledge, no comparable observation has been reported for 
eukaryotic Nav channels. NaChBac's preferential selectivity for sodium is reduced either by partial titration of its 
highly charged selectivity filter, when extracellular pH is lowered from 7.4 to 5.8, or by perturbation — likely ste- 
ric — associated with a nominally electro-neutral substitution in the selectivity filter (E191D). Although no single 
molecular feature or energetic parameter appears to dominate, our atomistic simulations, based on the published 
NavAb crystal structure, revealed factors that may contribute to the normally observed selectivity for Na over K. 
These include: (a) a thermodynamic penalty to exchange one K + for one Na + in the wild-type (WT) channel, in- 
creasing the relative likelihood of Na + occupying the binding site; (b) a small tendency toward weaker ion binding 
to the selectivity filter in Na-K mixtures, consistent with the higher conductance observed with both sodium and 
potassium present; and (c) integrated 1-D potentials of mean force for sodium or potassium movement that show 
less separation for the less selective E/D mutant than for WT. Overall, tight binding of a single favored ion to the 
selectivity filter, together with crucial inter-ion interactions within the pore, suggests that prokaryotic Nav channels 
use a selective strategy more akin to those of eukaryotic calcium and potassium channels than that of eukaryotic 
Nav channels. 



INTRODUCTION 

The classic experiments of Hille defined selectivity 
profiles for voltage-gated sodium (Nav) and potassium 
channels (Hille, 1971, 1972, 1973) , and later for nAChR 
channels at the neuromuscular junction (Adams et al., 
1980). These determinations were based on relative per- 
meability values calculated from shifts in reversal po- 
tential determined under voltage clamp. The approach 
offered a well-defined protocol for comparison of the 
ease with which different ions entered a particular chan- 
nel. Such measurements, systematically performed for 
alkali cations, and a series of variously sized organic cat- 
ions roughly defined the limits for the size, and to some 
degree the cross-sectional shape, of the narrowest part of 
each channel, henceforth known as the "selectivity filter." 
A similar analysis was performed for voltage-gated Ca 
channels from skeletal muscle (McCleskey and Aimers, 
1985) . Later, single-channel recording allowed direct 
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comparisons of unitary conductance, measured in the 
presence of different permeating ions, which provides 
an alternate indication of channel selectivity (Meech 
and Mackie, 1993; Thompson et al., 2009). Although, 
depending on the molecular details of conduction, rel- 
ative permeabilities and conductances obtained from 
these two approaches may differ, the analysis of reversal 
potential shifts resulting from external cation replace- 
ment endures as one convenient and widely used ap- 
proach to measure ion channel selectivity. 

Studies on Nav channels from frog node of Ranvier sug- 
gested a selectivity filter of asymmetric cross section, mea- 
suring a minimum of °^3 x 5 A (Hille, 1971 ) . Subsequently, 
the selectivity filter was later identified as the highly con- 
served signature motif, DEKA (see rSkMl pore sequence 
fragments in Fig. 1) , to which one amino acid residue was 
contributed by each of the homologous repeat domains 
of the Nav channel a subunit (Heinemann et al., 1992) . 
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In contrast, Hille's analysis for delayed rectifier K chan- 
nels suggested a symmetric selectivity filter of °^3 A in 
diameter (Hille, 1973) . The first crystal structure deter- 
mined for a K channel revealed a size closely matching 
this estimate, and showed that the lumen of the selectiv- 
ity filter was lined by backbone carbonyls (Doyle et al., 
1998) of residues from the K channel signature sequence 
(Heginbotham etal., 1994). In contrast, functional data 
for Nav and voltage-gated calcium (Cav) channels sug- 
gested a filter lumen lined by amino acid side chains 
providing a net negative charge (Kim et al., 1993; Tang 
etal., 1993; Yang etal., 1993). 

Functional characterization by the Clapham labora- 
tory of a homotetrameric, prokaryotic Nav channel, 
NaChBac (Ren et al., 2001), prompted hopes for an 
Nav channel crystal structure. After 10 years, these hopes 
were realized by Catterall and collaborators (Payandeh 
et al., 2011, 2012), using the NaChBac relative NavAb, 
closely followed by Zhang et al. (2012) and McCusker 
et al. (2012). Studies on these prokaryotic channels 
shattered the maxim that an asymmetric selectivity 
motif was required to achieve selectivity for Na + over K + 
and common divalent metal ions (Backx et al., 1992; 
Heinemann et al., 1992). In the prokaryotic Nav chan- 
nels, sequence alignment and structure are consistent 
with a selectivity ring made up of one glutamate residue 
from each of the four monomers, reminiscent of the 
EEEE selectivity filter of four-domain eukaryotic Cav 
channels (Fig. 1 D). In eukaryotic Nav channels, there 
is a conserved outer ring of acidic residues, three to 
four positions C terminal from the DEKA motif (Fig. 1 D) , 
but this outer ring appears to influence conductance 
and proton block, rather than determine selectivity (Khan 
et al., 2002). Although follow-up studies of NaChBac 
showed that it could easily be converted into a Ca- 
selective channel by the addition of one or more ad- 
ditional acidic residues into the outer vestibule-lining 
P loop (Yue et al., 2002), the detailed mechanism by 
which the NaChBac family of prokaryotic channels 
achieves Na selectivity with an EEEE selectivity filter has 
yet to be clarified. 

Similar to voltage-gated potassium channels, the 
NavAb channel is a homotetramer. Each monomer is 
composed of six transmembrane helices, S1-S6. The 
first four helices, S1-S4, form the voltage-sensing do- 
main, whereas helices S5 and S6, plus their connecting 
loop, form the pore domain, through which ions per- 
meate. Notably, the presumed selectivity filter of the 
NavAb channel is wider than those of K + channels. The 
resolution of the crystal structures renders unambigu- 
ous positioning of permeant ions difficult. Nevertheless, 
up to three binding sites were proposed based on a com- 
bination of functional and structural studies (Payandeh 
et al., 2011). T175 of each monomer forms the inner 
(cytoplasmic) site, S 1N , whereas main-chain oxygens of 
LI 76 aided by the water molecules comprise the central 



site, Scen, and the acidic side chain of E177 ligates perme- 
ant ion in the outermost, high field-strength site, Shfs 
(Fig. 1 A) . The question of how this filter gives rise to 
selective permeation, including the functional role of 
the varying polarity among its different binding sites, 
has yet to be resolved. 

The general physical principles underlying selectivity in 
prokaryotic and eukaryotic Nav channels may be simi- 
lar, but the molecular details show some striking differ- 
ences. A broad structural similarity is seen in the asymmetry 
along the pore axis. Both prokaryotic and eukaryotic 
Nav channels possess a selective periplasmic, or extra- 
cellular, section of the pore, with an outer vestibule that 
quickly grades into the narrow selectivity filter, followed 
by an expansive aqueous inner cavity, roughly aligned 
with the hydrophobic center of the surrounding bilayer, 
and then, a nonselective cytoplasmic mouth. The points 
of difference are tantalizing and include, in prokary- 
otes, the more negative net charge of the EEEE selectiv- 
ity ring and the radial quasi symmetry, which contrasts 
with the strongly asymmetric DEKA ring of eukaryotic 
Nav channels (Fig. 1 D). 

Among the eukaryotic members of the P-loop ion 
channel molecular family, Cav and Kv channels show 
functional characteristics that are most easily explained 
if each channel normally accommodates two to four 
conducting ions. These properties include tight bind- 
ing of the preferred ion species (Hess and Tsien, 1984; 
Morais-Cabral et al., 2001), strong coupling of unidi- 
rectional ion fluxes through the channel (Hodgkin 
and Keynes, 1955; Busath and Begenisich, 1982), and 
a nonlinear dependence of conductance on the mole 
fraction for mixed-ion solutions (Hille and Schwarz, 
1978). In contrast, eukaryotic Nav channels show a 
saturating dependence of single-channel conductance 
on sodium concentration with half-maximal conductance 
in the millimolar range (Hille, 1975b; Moczydlowski 
et al., 1984; French et al., 1994), very weak flux cou- 
pling (Busath and Begenisich, 1982), and other fea- 
tures consistent with only a single ion being bound to 
the channel most of the time (Ravindran et al., 1992; 
French et al., 1994) . Perhaps the most obvious connec- 
tion is that, functionally, bacterial Nav channels tend to 
be leakier to calcium ions than their eukaryotic coun- 
terparts, and several point mutations allow them to 
preferentially conduct calcium (Yue et al., 2002) . A key 
question that we address here is how bacterial Nav chan- 
nels can show selectivity and conductance closely ap- 
proximating those of their eukaryotic counterparts, 
whereas their multi-ion occupancy, suggested by avail- 
able structures, is more akin to that of eukaryotic Cav 
and Kv channels. 

We experimentally evaluate NaChBac's alkali cation 
selectivity under a variety of conditions and explore 
possible mechanistic and molecular underpinnings of 
that selectivity using computer simulations based on the 
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related NavAb crystal structure. Although these are two 
different channels, they display a considerable struc- 
tural similarity; for example, both channels are homo- 
tetramers with high sequence identity, especially in the 
pore domain. Judging by multiple studies of 2TM and 
6TM potassium channels, a wide variety of channels 
with apparent difference in sequence produces similar 
structures of the pore domain (within 1-2-A resolution 
for backbone atoms) . 

We identify possible determinants of conduction and 
alkali cation selectivity of bacterial sodium channels 
by analysis of molecular dynamics (MD) simulations 
and calculations of 1- and 2-D potentials of mean force 
(PMFs) for ions moving within the conducting pore. We 
find that the functions of selectivity and conduction re- 
semble those of eukaryotic Navs, even though patterns 
of ion interaction and occupancy are more like those of 
Cav and, to a lesser extent, Kv channels. 

MATERIALS AND METHODS 

Mutagenesis 

The original NaChBac channel construct in pTracer-CMV2 (Invi- 
trogen) was provided by D. Clapham (Howard Hughes Medical 
Institute, Children's Hospital, and Harvard University, Boston, MA). 
Single amino acid mutants were generated as described previously 
(Light et al., 2000) , using overlapping PGR amplification with oli- 
gonucleotide containing the sequence for the desired amino acid 
substitutions, followed by subcloning into pTracer (Invitrogen). 
All clones were completely sequenced. 

Electrophysiology 

Mammalian tSA 201 cells (Margolskee et al., 1993) were transfected 
with DNA encoding WT or mutant NaChBac, plus GFP using Poly- 
fect (QIAGEN). Significant channel expression occurred within 
24 h, after which whole-cell patch-clamp recordings were made at 
room temperature with an amplifier (Axopatch 200B; Molecular 
Devices). Patch pipettes were pulled from glass (8161; Corning; 
Potash-Rubium-Lead; softening temperature, 600°C; dielectric 
constant, 8.3) to a resistance of 1.5-2.5 MO.. Recordings were 
made 24 h after transfection in control external solution that con- 
tained (mM): 142.5 NaCl, 2 CaCl 2 , 2 MgCL, 10 glucose, and 10 
HEPES, pH 7.4. External ion replacement solutions were made by 
substituting 142.5 mM Na by Li, K, Rb, or Cs. Control intracellular 
(pipette) solution ("Na in ") contained (mM): 105 CsF, 35 NaCl, 
10 EGTA, and 10 HEPES, pH 7.2. "K^" was made by substituting 
35 mM NaCl by 35 mM KC1. For bi-ionic experimental conditions, 
the intracellular solution contained either 140 mM Na, K, or Cs. 
External solution changes were achieved by local superfusion of 
the replacement solution over the cell, with appropriate correc- 
tions for changes injunction potential (see Data analysis below). 

Data analysis 

Data were analyzed using Clampfit (Molecular Devices) and Igor 
(WaveMetrics) software. Peak I-V curves were fitted using I(V) = 
(V - V,„) * G max /(1 + exp((V half - V)/V slope )), where I is the mac- 
roscopic current, V is the command potential, V lev is the reversal 
potential, G max is the maximal conductance, V h:ll f is the half-activation 
potential, and V s i ope is the slope factor (mV/e-fold). For weakly 
permeant external ions, the reversal potential is expected to 
occur in a negative voltage range, at which very little conductance 
is activated (e.g., Fig. 4 C, left, internal sodium, external potassium). 



Thus, in the presence of weakly permeant ions, a prepulse to — 10 mV 
was applied to induce maximal activation, and V rev was determined 
from the instantaneous I-V relation, measured from the initial 
point in the tail current decay after steps to a series of voltages 
encompassing V lev . For weakly permeant ions, this protocol pro- 
vided a more precise and reproducible estimate of V rev . Relative 
permeabilities were calculated according to: Px/PNa = exp((E x — 
En 2 ) / (RT/F) ) , where Px is permeability to ion X, P Na is Na perme- 
ability, and RT/F is 25.4 mV. 

Net junction potentials were balanced to reduce the pipette 
current to zero before seal formation. In experiments where solu- 
tion replacement happened after the seal was established, the theo- 
retical junction potential for each solution pair was determined 
using JPCalc from Clampfit, which was also used for V rev correc- 
tion before the calculation of the permeability ratio. Series resistance 
compensation was applied conservatively to favor voltage-clamp 
stability; we note that the series resistance correction approaches 
zero near the reversal potential, the most critical measurement in 
this study. Net junction potential corrections applied to estimates 
of V rev fell in the range of 1.8-9.9 mV. 

All summary data are presented as mean ± SEM ( n) , where n is 
the number of determinations. Statistical significance was evalu- 
ated using the unpaired Student's t test; the criterion for a signifi- 
cant difference was taken to be P < 0.05, unless otherwise stated. 

Computational strategies 

The initial structure of NavAb was taken from the high resolution 
x-ray crystallographic structure with the Protein Data Bank acces- 
sion number 3RVY (Payandeh et al., 2011). The structure corre- 
sponds to a closed-pore conformation. The molecular simulations 
reported in this study were focused on the WT NavAb (Fig. 1 A) 
and the E177D mutant (Fig. 1 B). The tetrameric channel was em- 
bedded into a pre-equilibrated DMPC lipid bilayer and solvated in 
a 104 x 104 x 80-A 3 box filled with TIP3P water using CHARMM- 
GUI membrane builder protocol (Jo et al., 2008). The whole as- 
sembly was bathed with 150 mM NaCl, ensuring electroneutrality. 
Although Na + ions were not directly observed in the crystallo- 
graphic structure, three Na + ions were positioned in the filter as the 
starting conformation (as shown in Fig. 1 C), based on continuum 
electrostatic computations. All MD simulations were performed 
by the program CHARMM (Brooks et al., 2009) . The CHARMM-27 
force field was used for protein and lipids (MacKerell et al., 1998; 
Feller and MacKerrell, 2000), and Na + and K + ions as reported 
previously (Noskov et al., 2004; Noskov and Roux, 2008). 

The NpT ensemble was used for all simulations, with pressure 
set to 1 atm and temperature set to 315 K, Long-range electro- 
static interactions were treated by the particle mesh Ewald algo- 
rithm (Essmann et al., 1995). Nonbonded interactions were 
switched off at 10-12 A. The systems were simulated with periodic 
orthorhombic boundary conditions applied in all directions with 
the time step of 1 fs. After a staged equilibration with a gradual 
decrease in harmonic constraints that act on heavy protein atoms 
only, further nonconstrained equilibration was run for 23 ns. 

To unravel energetics of ion permeation, we used multidimen- 
sional umbrella sampling methods, a powerful computational 
technique used with considerable success in studies of K channels. 
Umbrella sampling simulations were performed with harmonic 
biasing potentials with a force constant of 10 kcal/ (mol ■ A 2 ) 
along the z axis. The zero position along the z axis is the center of 
mass of the selectivity filter backbone atoms of residues T175, 
L176, and E177. 1- and 2-D profiles for ions were computed. The 
lateral displacement of ion(s) was restrained to be within a cylin- 
der with a radius of 10 A and the central axis along the z axis. The 
reaction coordinate for each ion was the distance along the z axis 
between the ion and the center of mass of the selectivity filter of 
the protein, as just defined. The final snapshots of the conven- 
tional MD were used as the starting conformations for the umbrella 
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sampling. The sampling windows were spaced every 0.5 A from 
+14.5 to —10 A, resulting in 50 windows for 1-D PMF computations 
and 1 ,045 windows for 2-D PM Fs. The simulation time per window 
was set to 3 ns for 1-D PMF computations and 0.5 ns for 2-D PMF 
computations, respectively. 

1-D PMFs (e.g., Figs. 3 B and 6 B) were obtained by integration 
along a pathway for a single ion across a 2-D map (e.g., Fig. 3 C). 
Such 1-D PMFs, in general, would be expected to trace out a likely 
low energy pathway of one ion across the 2-D PMF surface, as illus- 
trated in Fig. 3 C, and would reflect the ion's interactions with its 
environment, including other ions, protein, and solvent. Conver- 
gence of computed 1-D PMFs obtained from a 2-D PMF map ob- 
tained after simulation times of 100-500 ps is illustrated in Fig. S3. 
The energy surfaces were rebuilt with the weighted histogram 
analysis method (WHAM) (Kumar et al., 1992; Grossfield, 2012), 
and the tolerance for WHAM was set to 0.001. 

In some cases, 2-D PMFs were further analyzed by extracting 
1-D "cross sections," for which one ion is held constant, whereas the 
other is moved on a single reaction coordinate across the profile 
(e.g., Fig. 3 E, and see Results). An indication of the significance 
of differences between PMFs or profiles of other parameters is 
given by SEMs estimated from standard block-average calculations 
(Rapaport, 2004) using progressively increasing block sizes (see 
Figs. 3, A and D, and 6, A, B, and D and associated text) . 

The single equilibrium dissociation constant ^(single) from 1-D 
PMF in the presence of a cylindrical constraint can be expressed 
as follows (Allen etal, 2004; Kim and Allen, 2011): 



where R is the radius of the cylindrical restraint oriented normal 
to the z axis with z„„„ = — 10 A and z, mx = 14.5 A. The w(z) was offset 
to zero for an ion in the bulk phase. 

The equilibrium dissociation constant for the double occupancy 
state of the channel, K D (double), can be expressed in accordance 
with Allen etal. (2004): 

Kp 1 (double) = jr , (2) 
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with the following integration limits: — 7 A < Z] < 20 A and — 10 A 

< z, < 20 A. 

Online supplemental material 

The supplemental figures provide additional experimental data 
to complement that presented in the main text, as follows. Fig. SI 
illustrates the time courses of wash-in and washout, as K + replaces 
Na + in the external solution. Combined with the results in the 
main text, Fig. 2 (A-C), the data make two points: (1) washout of 
K + , and its replacement by Na + , is significantly slower than the re- 
verse, suggesting stronger binding of K + under these conditions; 
and (2) washout of K + , in particular, depends significantly on the 
complement of intracellular ions, with the slowest being observed 
with intracellular solutions containing Cs + or Na + as the sole inter- 
nal alkali cation species. Fig. SI shows that no anomalous mole- 
fraction dependence is seen when external Na + is replaced by K + , 
with intracellular media containing 105 mM Cs + with 35 mM K + or 
Na + . This contrasts sharply with the observations for intracellular 
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Figure 1 . Pore domains 
for prokaryotic and eukary- 
otic Nav, Cav, and voltage- 
gated potassium channels. 

(A) Proposed locations of 
ion-binding sites in the se- 
lectivity filter of WT NavAb 
(Payandeh et al., 2011) and 

(B) the corresponding or- 
ganization of the selectivity 
filter in the NavAb E177D 
mutant. (C) NavAb crystal 
structure modeled within a 
DMPC lipid bilayer (see Ma- 
terials and methods) show- 
ing the permeation pathway 
lined by the selectivity filter. 
The nearest voltage-sensing 
domain and pore domain, and 
the farthest voltage-sensing 
domain, are removed for clar- 
ity. Chloride and sodium 
ions are shown as cyan and 
yellow spheres, respectively. 
The lipid phosphate groups 
are represented by orange 
spheres. The simulated aque- 
ous phase is bounded by the 
gray lines. (D) Comparison 
of pore domain sequences 
among prokaryotic and eu- 

karyotic voltage-gated cation channels (PI, pore helix 1; SF, selectivity filter; P2, pore helix 2). Ident, percent identity; Simil, percent 
aligned similar residues compared to NaChBac. Note that Qjf, the net charge on the selectivity filter, includes both inner and outer 
charged rings in the eukaryotic sodium and calcium channels. The rat skeletal muscle sodium channel, rSkMl, is also known as rNavl.4. 



NaChBac 


173 


NavAb 




159 


NavRh 




161 


rSkMl- 


-DI 


382 


rSkMl- 


-DII 


737 


rSkMl- 


-Dili 


1119 


rSkMl- 


-DIV 


1510 


hCavl 


2-DI 


345 


hCavl 


2-DII 


688 


hCavl 


2-DIII 


1117 


hCavl 


2-DIV 


1446 



Shaker 
KvAP 



423 
177 



YFGNLQLSLLTLFQVVTLESWASGVMRPIFAEVWWS 
WFGTLGESFYTLFQVMTLESWSMGIVRPLMEVYPYA 
SRWDLGISLITLFQVLTLSSWETVMLPMQEI YWWSW 

SYDTFSWAFLALFRLMTQDYWENLFQLTLRAAGKTY 
HMNDFFHSFLIVFRILTGEWIETMWDCMEVAGQAMC 
NYDNVGLGYLSLLQVATFKGWMDIMYAAVDSREKEQ 
FNFETFGNSICLFEITTSAGWDGLLNPILNSGPPDC 

NFDNFAFAMLTVFQCITMEGWTDVLYWMQDAMGYEL 
TFDNFPQSLLTVFQILTGEDWNSVMYDGIMAYGGPS 
DFDNVLAAMMALFTVSTFEGWPELLYRSIDSHTEDK 
NFQTFPQAVLLPLRCATGEAWQDIMLACMPGKKCAP 

NSFFKSIPDAFWWAVVTMTTVG3TGDMTPVGVWGKIV 
NSSIKSVFDALWWAVVTATTVGYGDVVPATPIGKVI 



208 


100 


100 


-4 


194 


61 


79 


-4 


196 


55 


73 


-4 


417 


31 


69 


-5 


772 


0 


0 




1154 


0 


0 




1545 


0 


0 




380 


41 


64 


-8 


723 


47 


60 




1152 


36 


54 




1478 


0 


0 




458 


40 


60 


0 


212 


0 


0 


0 



160 Prokaryotic Na channels: Alkali cation selectivity 



solutions containing either only K + or only Na + , which reveal a 
dramatic anomalous mole-fraction dependence (Fig. 4 D) . Fig. S2 
provides voltage-clamp records from NaChBac E191D taken at 
extracellular pH values of 7.4 and 5.8. Even though the E191D 
mutant is much less selective for Na + over K + than is the WT channel, 
an additional significant loss of selectivity (i.e., an increase in 
Pk/Pn's) was observed with a decrease of extracellular pH (Fig. 5 B) . 
Fig. S3 illustrates convergence of a 1-D PMF computed from a 2-D 
map. The online supplemental material is available at http:// 
www.jgp.org/cgi/content/full/jgp.201311037/DCl. 

RESULTS 

Overview 

In this Results section, we stress the overriding goal of 
our study: to systematically integrate experimental and 
computational results that suggest the mechanisms un- 
derlying selectivity among alkali cations by bacterial Nav 
channels, exemplified by NaChBac and NavAb. To this 
end, we interleave results from experiments and simula- 
tions throughout. We begin with the classic approach of 
measuring relative permeabilities, based on reversal po- 
tential shifts observed when the external cation species 



is changed, and calculations of the PMF as a single ion 
is moved through the channel. Intuitively, one might 
expect that reversal potential shifts generated by exter- 
nal ion substitution are expected to be primarily influ- 
enced by the energetics of ion entry into the pore (Patlak, 
1960; Bezanilla and Armstrong, 1972), and hence by 
the outermost parts of the pore structure. It is this part 
of the pore structure that is likely to differ least between 
open and closed conformations in channels for which 
control of channel opening is presumed to lie at the S6 
bundle crossing on the cytoplasmic side of the selectiv- 
ity filter. Thus, we believe that the putative preopen 
structure of Payandeh et al. (2011) is an appropriate 
basis for a consideration of selectivity based on relative 
permeability measurements. 

Subsequently, we present experimental results observed 
after changes in the ion species on both extracellular and 
intracellular sides of the membrane. These experiments 
are more likely to yield results indicating whether or not 
ion-ion interactions play a prominent role in ion selec- 
tivity and conduction. Their analysis required simulations 
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Figure 2. NaChBac relative permeabilities 
to alkali cations show the same sequence as 
for eukaryofic Nav channels, with small quan- 
titative differences for some individual ions. 

(A) Whole-cell NaChBac-mediated currents 
in the presence of different extracellular 
alkali cations. Current amplitudes for the 
substituted ions are normalized to the maxi- 
mum inward sodium current in the same cell. 

(B) Representative I-V relationships for each 
cation substitution shown in A. Ordinate val- 
ues show current amplitudes with 142.5 mM 
of external Na; for other ions, currents are 
normalized as in A. (C) Relative permeabili- 
ties, from I-V reversal potential shifts after ex- 
ternal cation substitutions, calculated using 
the equation shown. Closed bars, with error 
bars, correspond to NaChBac data, and open 
bars show relative permeability of frog node of 
Ranvier Na channels, as reported by Hille (1972) . 
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in which at least two ions occupy the channel, which led 
us to explore the implications of the resulting 2-D PMFs 
that emerge (see below) . 

Relative permeability fingerprint 

The qualitative selectivity fingerprint (Fig. 2), defined 
by the sequence of relative permeability values, P x , where 
P Na is defined as 1, matches both the sequences deter- 
mined from similar experiments for eukaryotic Nav chan- 
nels (e.g.,Hille, 1972) and the sequence for equilibrium 
binding to high field-strength anionic sites (Eisenman, 
1962). Our experiments used 105 mM Cs as the major 
internal cation with 35 mM of internal Na to allow easy 
measurement of either positive or negative shifts in Erev, 



with a quasi-physiological inward electrochemical gradi- 
ent for Na over most of the working experimental volt- 
age range. The resulting selectivity sequence is P L i ~ Pn 3 
> Pr ~ PRb ~ Pes- 

Moving ions through the pore: PMFs for single 
and double occupancy 

As a first step toward defining the energetics of ion per- 
meation, we evaluated TD PMFs as single ions were 
moved through the NavAb model pore (Fig. 3 A) . In the 
absence of other ions, a solitary Na + would bind much 
too tightly to permit physiological rates of throughput 
(single-channel conductance, ^T2 pS or ^lO 7 ions/s~\ 
with a 100-mV driving force; Ren etal., 2001). From the 




-10 -5 0 5 10 15 



Reaction Coordinate (A) 




-5 0 5 10 15 20 




-5 0 5 10 15 20 

Reaction Coordinate (A) 



B 




-25 J , , , , , r 

-10 -5 0 5 10 15 

Reaction Coordinate (A) 




Reaction Coordinate (K + ) (A) 

4 -i 




■•• K-Na 



0 5 10 15 20 

Reaction Coordinate (A) 



Figure 3. Calculated energy profiles and 
contour maps, for single-ion and two-ion 
pore occupancy, which determine the 
movement of ions across the NavAb pore. 
(A) 1-D potential of PMF for the move- 
ment of a single Na + ion across the selec- 
tivity filter of NavAb. Note that single-ion 
binding would be much too tight to per- 
mit the experimentally observed conduc- 
tances. (B) Superposition of 1-D PMFs for 
movement of a single Na + ion in the cases 
of single and double occupancy (see low 
energy pathway, i-v in C, Na-Na case.). 
For double occupancy, the barrier for 
ion exit back to the extracellular solution 
is drastically reduced, reflecting weaker 
binding and allowing larger conductance. 
(C) 2-D PMFs for the movement of two 
ions, Na + and/ or K + , as indicated, with la- 
beled points ( i-v) shown along the lowest 
energy permeation pathway. Na-Na and 
K-Na describe the plots where a Na + , or 
a K + ion, moves in the presence of a Na + 
ion. K-K and Na-K are plots where K + or 
Na + moves in the presence of a K + ion. In 
every case, both ions were free to move, as 
indicated in the insets. (D) 1-D PMFs for 
different double-occupancy combinations 
from the 2-D contour plots shown in C: 
Na-Na, K-K, K-Na, and Na-K. Each pro- 
file was obtained by integration from the 
2-D PMF maps (see Materials and meth- 
ods). (E) 1-D cross sections of 2-D PMFs 
for different double-ion occupancy of the 
selectivity filter, when one ion is fixed at 
the position of binding site ii, (see C). The 
cross sections illustrate the energy scape 
experienced by the second ion-entering 
selectivity filter. As is apparent from the 
differences between D and E, significant 
energetic differences arise depending on 
the particular placement, order of entry, 
and pairing of the different ions. 
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TABLE 1 



Gibbs free energies and equilibrium dissociation constants for single- and double 


-ion occupancy from 


2-D PMF computations 


Occupancy group 


Occupancy 


Dissociation constant (mol/L) 


AG (kcal/mol) 






K D (single) 


K D (double) 


AG 


Single-ion occupancy 


Na 


4.88E-14 




-19.19 




K 


7.56E-12 




-16.03 


Na at innermost position 


Na-Na 




4.27E-03 


-3.42 




K-Na 




8.44E-03 


-2.99 


K at innermost position 


K-K 




9.85E-03 


-2.89 




Na-K 




7.29E-03 


-3.08 



Values were calculated from the following data: Ko (single) for Na, from Fig. 3 A, and from equivalent data for K K D (double) from the data in the 2-D 
PMFs in Fig. 3 G. 



one-ion PMF, dissociation rates of 10 1 ~ 10~ 2 ions/s 1 
are predicted from a single-barrier calculation based on 
Eyring rate theory. This is reminiscent of tight binding 
of Na + and Ca 2+ to binding sites in model channels 
(Nonner and Eisenberg, 1998; Krauss et al., 2011). To 
reconcile large ionic fluxes, Krauss et al. (201 1 ) argued for 
the existence of multiple binding sites and the impor- 
tance of the local crowding effects. More realistic perme- 
ation pathways appear on the energy landscape when 
two monovalent cations enter into the pore (Fig. 3 C) , 
with energy minima for four doubly occupied configu- 
rations in the range of —2 to —4 kcal/mole (Fig. 3, B 
and D). Then, low energy permeation pathways (blue 
tones) become clearly apparent; the potential profiles 
traversed by a single ion were obtained by integrating 
along the reaction coordinates (Fig. 3 D and see Materi- 
als and methods, Computational strategies) . Similar find- 
ings have been reported recentiy from computational 
studies by Furini and Domene (2012) and by Corry and 
Thomas (2012) . They also emphasize the need for a sec- 
ond bound ion to achieve physiological conductance 
levels. The relatively deep energy well for Na-Na would 
contribute to higher stability of Na-Na-occupied filter, 
with the essentially flat entrance barrier for Na-Na also 
playing a significant role (Z = 10 ^ 12 A) in discrimina- 
tion among different ions. It is important to stress that 
uncertainties in 2-D PMF computations are rather high, 
but the resolution in the computed maps is still suffi- 
cient to clearly discriminate between Na-Na and K-K 
occupancy in the filter. Interestingly, in the case of K-K 
binding to the filter, there is a small barrier separating 
two ion-binding energy minima. 

Extended analysis of the 2-D PMFs revealed multiple 
binding sites separated by small barriers (indicated in 
Fig. 3, C-E) . The integrated single-ion PMFs (e.g., Fig. 3 D) 
and TD "cross-section" profiles (e.g., Fig. 3 E), derived 
from the 2-D PMF surfaces, help us to understand the 
significance of this apparent complexity in the filter. We 
estimated the corresponding equilibrium dissociation 
constants in the presence of constraints for the singly 
and doubly occupied channel using the following theo- 
retical developments of Kim and Allen (2011) andDerebe 



et al. (201 1 ) . In general, the channel binds Na + in prefer- 
ence to K + , but when two ions are bound, the magnitude 
of Na/K selectivity for binding at the first site depends 
on the identity of the second ion. With symmetric oc- 
cupancy of the filter, binding of an Na + -Na + pair is fa- 
vored over a K + -K + pair (Table 1 ) . Up to three well-defined 
potential binding sites can be distinguished, further sup- 
porting the idea that multiple ions can occupy the pore. 
The free energy cost to bring in a third ion would render 
this occupancy unstable, thus leading to a permeation 
event. The effective free energy cost caused by ion-ion 
repulsion for bringing a second ion as determined from 
the single- and double-occupancy dissociation constants, 
Kg and Kj, is ^13-16 kcal/mol depending on the cat- 
ion species. 

The mixed occupancies, Na-K and K-Na, yield inter- 
mediate depths, reflecting an effect on the stability of 
the binding around the middle of the selectivity filter. 
However, the PMF calculations are accurate only within 
^1 kcal/mol according to our block-average analysis plac- 
ing restraints on the interpretation of obtained PMFs. 
Nevertheless, in the limits of method resolution, the 
barriers at the extracellular side of the channel display 
dependence on the combination of ions present, al- 
though the differences are subtie. In the case of the bar- 
riers, the mixed ion combinations provide the extremes, 
whereas double occupancy, by a single-ion species, Na- 
Na and K-K yields the extremes. The number and the 
conformation of binding sites definitely differ consider- 
ably depending on the identity of the second ion. The 
free energy profile for Na + , double occupancy (the solid 
black line in Fig. 3 D) , has a local minimum of around 
Z = 3 A, which is not present in the case of K + (the solid 
red line in Fig. 3 D). Also, binding sites around Z = 4.5 A 
for Na + shifted to the right by M).5 A with respect to those 
for K + . 

"Selectivity" is asymmetric and shows a complex 
dependence on ionic composition of the solutions 

Garber (1988) reported a striking asymmetry of reversal 
potentials for skeletal muscle Nav channels in planar lipid 
bilayers, measured when a bi-ionic Na/K gradient was 
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reversed. Similar behavior has also been reported re- 
cently by a mutant version of nonselective NaK channel 
byDerebe et al. (2011). 

We observed an analogous result for NaChBac whole- 
cell currents, when K + ext was exchanged for Na + ext in the 
presence of either internal Na + or K + (Fig. 4 C) . The mag- 
nitude of the reversal potential shift, after external ion 
replacement, was ^100 mV with internal Na + but only 
^40 mV with internal K + , yielding a 10-fold difference 
in Pk/Pnb calculated with oppositely directed gradients 
(Pic/Pi^ = 0.2 with internal K + , and 0.02 with internal 
Na + ). In bi-ionic conditions, with isotonic Na + and K + 
solutions on opposite sides of the membrane, the rever- 
sal potentials were approximately +45 mV (Na + inside) 
and —100 mV (Na + out; instantaneous I-V relations; data 
points plotted as open triangles; see Fig. 4 C) . This func- 
tional asymmetry may reflect a common property among 
the various members of the P-loop family of channels, 
which possess a nonselective entrance to the intracellular 



cavity, and a selectivity filter located nearer to the extra- 
cellular mouth. 

In doing these experiments, we also noted a kinetic 
asymmetry. The approach to a new steady-state level of 
peak currents elicited by a succession of identical voltage- 
clamp pulses took significantly longer for washout of ex- 
ternal K + by Na + than vice versa (Fig. 4, A and B) , suggesting 
stronger binding of K + within the pore. This difference 
was most obvious with Na + as the only internal cation spe- 
cies, but it was also significant with K + as the only intracel- 
lular cation. Slow washout of potassium was also seen with 
cesium as the only internal cation, but it was not detect- 
able using mixtures of cesium or potassium with sodium 
(Fig. SI, A-D). 

Dependence of conductance on extracellular mole fraction 
in sodium-potassium mixtures 

Progressive substitution of external Na + by K + reveals a 
nonmonotonic (anomalous) dependence of conductance 
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Figure 4. K + binding in the pore, asymmetry 
of bi-ionic Na-K reversal potentials, and anom- 
alous mole-fraction dependence of conduc- 
tance. (A) Representative wash-in and washout 
time courses, observed when a 142.5-mM Na + 
extracellular solution is fully replaced by 
142.5 mM K + , and is subsequently washed 
out. Insets show current traces in response to 
a — 10-mV stimulation every 10 s with [Na + ]j = 
140 mM. (B) Time constants derived from 
exponential fits to the wash-in (Na— >K) and 
washout (K— >Na) for four different intracel- 
lular solutions. Internal alkali cation concen- 
trations were (mM): 35 Na (35 Na and 105 
Cs), 35 K (35 K and 105 Cs), 140 Na (140 Na), 
and 140 K (140 K). (C) I-V relationships indi- 
cating bi-ionic reversal potentials determined 
for I^/Nai (left) and Na 0 /Kj (right). The re- 
versal potentials change asymmetrically when 
the ionic gradients are reversed in direction. 
Note that clear measurement of reversal po- 
tentials for the bi-ionic case internal sodium/ 
external potassium required instantaneous 
I-V data (left, open triangles; see Materials 
and methods, Data analysis, for more details). 
In all other cases, reversal potentials were de- 
termined from peak I-V relations. (D) Maxi- 
mal conductance, Gmax, normalized with 
respect to the control value in full sodium 
external solution (/ Na = 1), plotted against 
the extracellular mole fraction of Na (^ a ). 
Prominent maxima are seen at intermediate 
values of ^ a for both internal Na and internal 
K, and a small but significant local minimum 
(circled) appears, for internal Na only, near 
/ Xa = 1. Note that error bars (±SEM) are plot- 
ted on all points. See text for further details. 
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on mole fraction of sodium, (Fig. 4 D). Such behav- 
ior has frequently been associated with multi-ion occu- 
pancy of channels and repulsive interactions among ions 
in the channel, although alternate explanations have 
been proposed. The following features are of interest in 
our data. With a single monovalent cation species inside — 
either sodium or potassium — there is a general tendency 
for conductance to increase in a supra-linear fashion up 
to fx* = 0.6. With potassium as the sole internal cation 
species, a conductance maximum occurs at fa ~ 0.8, 
and the conductance then declines monotonically to the 
value characteristic of a full sodium external solution. 
Notably, the appearance of this maximum appears spe- 
cific to interactions between Na + and K + , as we saw no 
maximum with internal solutions containing 105 mM 
Cs + , plus 35 mM of either Na + or K + (see Fig. SI, E and F) . 
If the sole internal cation species is sodium, a shallow 
local minimum also appears in the conductance at^j a ~ 
0.95 (Fig. 4 D). This result provides further evidence 
that the last K + ion to exit binds tightly enough to make 
it relatively difficult to completely wash potassium ions 
out of the channel (see also Fig. 4 A) , consistent with 
the tight binding predicted by TD PMFs. 

Two features of these data deserve attention. First, the 
maxima we see for intermediate mole fractions (/n 3 = 
0.6-0.8) are unusual, and to our knowledge have not 
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been reported previously for other systems. Thus, the 
conductance is higher for external mixtures of Na and 
K than with either Na or K as the sole external species. 
This would be consistent with an evolutionary optimiza- 
tion for conduction in physiological solutions contain- 
ing mixtures of sodium and potassium. 

In other reported cases of an anomalous mole-fraction 
effect involving K or Cav channels, a minimum in mole- 
fraction dependence has been seen, suggesting that one 
ion species binds strongly enough to impede permeation 
by another (Hille and Schwarz, 1978; Hess and Tsien, 
1984; Eisenman etal., 1986). A similar mechanism may 
underlie the local minimum that we show in Fig. 4 D, 
at f N3 ~ 0.95, with Na + as the internal monovalent cat- 
ion species. 

Simulating ion-ion interactions 

To explore the role of a particular bound ion on the 
entrance barriers for the second ion, we used a cross- 
section methodology from Li et al. (2010) applied to a 
bi-ionic solution, where a TD PMF ("cross section") for 
a permeating ion, in the presence of a fixed coperme- 
ant ion, was obtained from a 2-D map (Fig. 3 E) . The 
bound ion was fixed at Z = 1.0 A. As expected, the barri- 
ers presented near the entrance (Z > 12 A) display rela- 
tively litde dependence on a particular combination of 



Figure 5. Roles of side-chain 
length and negative charge 
density in the selectivity filter. 
Both experimental manipula- 
tions (E191D substitution and 
lowering of the external pH) re- 
duce selectivity among alkali cat- 
ions. (A) Representative E191D 
mutant-mediated currents in the 
presence of different extracel- 
lular alkali cations at pH 7.4. 
The bar graph on the right shows 
significant increases in both 
Ps/PNa and P Rb /P Na after the 
E191D substitution. (B) NaCh- 
Bac WT current traces recorded 
in the presence of Na and K at 
pH 5.8. To better resolve rever- 
sal potentials for external K, 
they were determined from in- 
stantaneous I-V relations after 
a depolarizing prepulse used to 
maximally activate the channels 
(see Materials and methods, and 
Fig. 4). Permeability ratios for 
WT and E191D mutant channels 
at neutral and acid pH are sum- 
marized in the bar graph on the 
right. Discrimination between Na + 
and K + is reduced at pH 5.8 for 
both NaChBac WT and for the 
E191D mutant. Probabilities, P, 
for different comparisons are 
indicated as follows: *, <0.05; 
**, <0.01;*** <0.0001. 
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cations. Within the selectivity filter (Z = 4-8 A) , the situa- 
tion is drastically different among the Na-Na, Na-KK-Na, 
and K-K combinations. The interaction between cat- 
ions in the filter can significandy affect the relative height 
of the barrier and the depth of the well experienced by 
a second cation. The ion bound to the lowest affinity 
site (Z = 1.0 A) also affects the location of the binding 
site for the entering cation. Notably, K-Na mixtures 
tend to display flatter energy profiles, suggesting that 
the conductance for a mixture of cations could be as 
high as for pure Na + solution based on reduction of the 
entrance barrier for the second cation. 

It is important to stress some limitations of our pres- 
ent analysis. It is based on a crystal structure for closed 
state, where complete integration of the permeation 
profile is impossible. Estimates of the mutual repulsion 
among permeant ions, and corresponding pore-ion in- 
teractions, would also depend on the force field used to 
analyze them. 

Permeation volume of the selectivity filter 
and discrimination among alkali cations 

Discrimination among alkali cations by NaChBac is 
qualitatively similar to that shown by eukaryotic Nav chan- 
nels (Hille, 1972). For NaChBac, under the conditions 
of our experiments (single monovalent cation species, 
X, in the extracellular solution, with an intracellular so- 
lution containing 35 mM Na and 105 mM Cs) , the relative 
permeabilities, Px/Pns, fall into the following sequence: 
Li ~ Na > K ~ Rb ~ Cs. Given the difference in net charge 
on the selectivity rings between bacterial (EEEE) and 
eukaryotic sodium channels (DEKA), an obvious ques- 
tion arises: how do changes in the net charge on the 
selectivity filter affect selectivity? 

Replacement of the selectivity filter glutamate with 
aspartate (E191D in NaChBac and E177D in NavAb) 
would yield a DDDD selectivity ring. With minimal per- 
turbation in the backbone structure, this substitution 
might have been expected to lead to a small increase in 
diameter, with a consequent reduction in negative charge 
density at the selectivity filter, and because of the shorter 
aspartate side chains. Actually, the energy-minimized 
structure for E177D shows a slightly reduced diameter 
because of the more direct projection of the side chain 
toward the pore axis than for the longer, more flexible 
glutamate (Figs. 1 B and 6 A) . The resultant decrease in 
radius is associated with the removal of the multiple en- 
ergy minima seen in 2-D PMFs, and it yields almost indis- 
tinguishable single-ion PMFs for Na + and K + (Fig. 6 D) . 

Experimentally, this mutant showed a decrease in se- 
lectivity for Na + , with significant increases in relative 
permeability to both potassium and rubidium compared 
with the WT channel (Fig. 5 A) . In a complementary set 
of experiments, lowering the extracellular pH from 7.4 
to 5.8, which is expected to titrate some of the charge 
on the EEEE selectivity ring, reduced the discrimination 



among alkali cations (Figs. 5 B and S3). For both WT 
NaChBac and the mutant E191D, the lowered external 
pH yielded an increased relative permeability to potas- 
sium. The permeability increase was slightly larger for 
the E191D mutant, such that this construct did not distin- 
guish between extracellular Na + and K + at pH 5.8. Collec- 
tively, the experiments in this section suggest that charge 
density within the selectivity ring is an important con- 
tributor to overall selectivity among the alkali cations. 

To test whether the pore cross section is affected by 
E191D replacement, we ran a complementary set of all- 
atom MD simulations on a mutated pore. The calcu- 
lated cross-sectional profiles for NavAb E177D are shown 
in Fig. 6 A. The pore radius at the constriction is still 
considerably wider than that for K channels. This find- 
ing is in keeping with previous data on the ion hydra- 
tion in the pore (Corry and Thomas, 2012; Treptowand 
Klein, 2012), which suggest that the permeant cation 
(either Na + or K + ) is able to maintain partial hydration 
and could possibly be accommodated by either WT or 
E/D mutant. The pore radius decreased slightly with 
the E177D mutation despite the shorter side chain 
introduced. To further investigate the mechanisms by 
which the E/D mutation might reduce selectivity, we 
again used a 2-D PMF-based strategy similar to that in 
the WT sections. The resulting 2-D PMFs for Na + /Na + 
and K + /K + are shown in Fig. 6 C, and the corresponding 
single-ion PMFs are shown in Fig. 6 B. The PMF calcula- 
tions show considerable changes in the energy surface, 
which controls ion transport across NavAb E177D. The 
computations suggest that the degree of ion stabiliza- 
tion in two of the three proposed binding sites present 
in the WT filter (Fig. 1 A) is significandy affected by the 
mutation. The first binding site (S H fs) is destabilized for 
both Na + and K + (Fig. 6 B), whereas the second binding 
site (Si N ) remains virtually unaffected. It is tempting to 
think that E/D mutation not only affects the cross sec- 
tion of the pore but also ion coordination environment 
and number of sites available to an ion along the per- 
meation pathway. The significant destabilization of one 
of the ion-binding sites caused by mutation would re- 
duce the likelihood of two ions binding simultaneously 
to the filter. 

The computed PMFs suggest that the differences in 
Na/K profiles, present in WT, are largely gone when this 
mutation is present, and it is also evident that there are 
several significant rearrangements in the ion coordina- 
tion in the filter region. The entrance binding site, ii, 
presented in the WT (Fig. 3 C) is eliminated in the mutant 
for both Na + and K + binding to E/D-NavAb, whereas 
central binding sites (Hi and iv) display similar stability 
for Na + and K + . This finding is in agreement with our 
experimental data on ion selectivity (Fig. 5 A) . Thus, the 
reduced selectivity of the E-D mutant may be primarily 
attributed to the removal of one of the binding sites from 
the selectivity filter coupled with a nondiscriminatory 



166 Prokaryotic Na channels: Alkali cation selectivity 



binding of Na + and K + to the remaining site (Fig. 6, B 
and C) . These results highlight an importance of the 
multiple binding sites present in the selectivity filter, 
which appears to modulate selective conductance. 

DISCUSSION 

Overview: Defining selectivity 

A definition of the "selectivity" of an ion channel could 
be based on measurements of ion-binding affinity, con- 
ductance, or reversal potential for currents measured 
under particular conditions. Under fairly general condi- 
tions (Patlak, 1960), "relative permeabilities" thus calcu- 
lated (see Materials and methods) direcdy reflect the rate 
of unidirectional flux across the membrane (Goldman, 
1943; Hodgkin et al., 1952), or the rate of entry into the 
channel (Hille, 1975a). Although this simple physical in- 
terpretation likely breaks down with complex barrier pro- 
files or ion-ion interactions (Eisenman and Horn, 1983) , 
relative permeability provides a widely used method to 



establish selectivity sequences among the various ions 
that can pass through a channel. 

In the following discussion, we establish selectivity se- 
quences among the alkali cations for NaChBac and its 
point mutant, E191D, at both normal and acidic pH. 
We go on to examine the complex concentration depen- 
dence of conductance to gain insight into ion binding 
to the channel and intrapore ion-ion interactions among 
ions within the pore. In the context of these experimen- 
tal results, we use MD and free energy calculations based 
on the structure of NavAb, a close relative of NaChBac. 

Summary of experimental results 

Homotetrameric NaChBac, despite its more symmetric 
structure, qualitatively resembles eukaryotic Nav chan- 
nels in that it selectively conducts Na in preference to 
most other monovalent ions and common physiological 
divalents. The selectivity "fingerprint," defined as the 
order of relative permeabilities, Px/PNa, for the alkali 
cations is indistinguishable from that seen for eukaryotic 
channels: Li ~ Na > K ~ Rb ~ Cs (Fig. 2 C) . As suggested 
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Figure 6. Predicted consequences of a 
charge-conserving selectivity filter muta- 
tion (NavAb E177D). (A) Compared with 
NavAb WT, there is a decrease in radius 
of the cylindrical cross section available 
for ion permeation along the pore axis 
for the E177D mutant. (B) 1-D PMFs inte- 
grated from 2-D PMF maps (c.f., Fig. 3 D) 
from C for Na + and K + transport through 
NavAb E177D. The lower entry barrier for 
K + would favor a higher value of Pk/Pns, 
thus lowering selectivity, as experimentally 
observed for NaChBac E191D. (C) Con- 
tour plots of 2-D PMFs, which provide the 
1-D plots in D, for occupancies by Na + /Na + 
and K + /K + in the NavAb E177D selectivity 
filter. (D) Average coordination numbers 
for sodium or potassium ions at different 
ion positions in the transport process for 
NavAb WT and for E177D are shown as 
solid and dotted lines, respectively (for 
coordination group identity, see the color 
key in the figure) . Coordination numbers 
represent the mean numbers of oxygen 
atoms, which are found within 3.5 A of 
an ion, as water, carbonyl, carboxylate, 
or hydroxyl groups. The choice of the 
unique confinement radius is a challeng- 
ing issue. We based our choice for the size 
of the coordination shell on an analysis that 
was chosen in accordance with our previ- 
ous studies on ion hydration by aqueous 
clusters and model protein sites (Yu et al., 
2010; Lev et al., 2013), where it was found 
to faithfully reproduce relative free ener- 
gies. In addition, the total number of co- 
ordinating oxygen atoms is shown by the 
solid red and black lines. 
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by earlier work (Ren et al., 2001 ) , the critical physiologi- 
cal selection for Na over K appears to be somewhat 
more stringent for NaChBac than for eukaryotic chan- 
nels (NaChBac, Pk/P N; , = 0.03 + 0.007(4) compared with 
0.048-0.086 for four eukaryotic Nav channel variants; 
see Table 14.2 in Hille, 2001). 

A functional asymmetry, common to NaChBac and 
eukaryotic channels, is seen in the inequality of magni- 
tude of reversal potentials observed when a bi-ionic Na/K 
gradient is reversed. Reversal potentials of approxi- 
mately — 100 mV (Nai/Iv,) and approximately +40 mV 
(Ki/Na„) indicate a 10-fold stronger discrimination be- 
tween K and Na when the internal cation species is Na 
than for the reverse case (Fig. 4 C). Alkaloid-modified 
rat muscle Nav channels (bound by a large alkaloid, ba- 
trachotoxin or veratridine, in the cytoplasmic cavity) 
show many of the qualitative properties of unmodified 
channels, including an asymmetry of bi-ionic reversal 
potentials (Garber, 1988) . This suggests that important 
determinants of selectivity reside near the outer end of 
the pore. The results also indicate the usefulness of a 
"pre-open" structure as a relevant basis for simulations 
of selectivity, particularly discrimination among extra- 
cellular ions, even though the intracellular bundle- 
crossing gate may be closed or modified. 

The dependence on external mole fraction in bi-ionic 
Na-K mixtures (f N2l = [Na] 0 /([Na] 0 + [K] Q )) of the nor- 
malized maximal conductance is highly nonmonotonic 
(Fig. 4 D). A clear maximum appears at intermediate 
values of/v a (0.6-0.8), suggesting that mixed occupancy 
by Na and K, as would be expected under physiological 
conditions, enhances conductance. In addition, a less 
conspicuous but significant local minimum, which occurs 
at low levels of K,, (f Nli ~ 0.95-0.97), hints that a single 
potassium ion binding at the outer end of the pore may 
inhibit sodium permeation. Anomalous (nonlinear and 
nonmonotonic) mole-fraction dependence has generally 
been associated with multi-occupancy with ion-ion in- 
teractions in the pore. This is consistent with our own cal- 
culations (Fig. 3) and those of others (Yue and Marban, 
1990), but alternate explanations have been proposed 
(Nonner et al., 1998; Gillespie and Boda, 2008). NaCh- 
Bac 's strong anomalous mole-fraction dependence, with 
its distinctive enhancement of conductance at interme- 
diate mole fractions, was observed only with intracellu- 
lar solutions containing solely Na, or solely K, as cations. 

Two maneuvers, which might be expected to change 
the charge density on the selectivity filter (E191D muta- 
tion and lowering of pH), each reduce NaChBac 's abil- 
ity to discriminate between Na and K (Fig. 5) , as predicted 
by Corry and Thomas (2012). Accordingly, our data 
show that lowering pH from 7.4 to 5.8 increases Pk/Pnsi 
significantly from 0.03 + 0.006(4) to 0.15 ± 0.03(6) (P = 
0.018; see Fig. 5B). 

The E191D substitution also reduces selectivity, with 
Pr/Pi*. = 0.57 ± 0.02(4) at pH 7.4 and a further increase 



to 0.71 ± 0.01 (4) at pH 5.8. It is not intuitively obvious 
why the decrease in pore diameter (Fig. 6 A) should be 
associated with a decrease in selectivity, but at least it 
does suggest that energetic factors beyond simple steric 
filtering, or simple binding interactions dependent on 
the field strength near a coordinating site, must be con- 
sidered to account for this change, together with the 
other complexities of the dataset. 

Computational studies 

Insight from atomistic simulations is essential to study 
the mechanism of ion permeation through NavAb chan- 
nels. Without such a quantitative approach, coordinated 
with experimental studies, the binding sites proposed 
based on the crystal structure would remain speculative. 
The deep energy wells shown in 1-D PMF profiles (e.g., 
Fig. 3 A) indicate that single Na + or K + ions are extremely 
unlikely to permeate though the channel. Calculations 
based on 2-D PMFs provide more realistic permeation 
pathways, with energy minima in the range of —2 to —4 
kcal/mol (Fig. 3, B and D) . The blue tones in the en- 
ergy landscapes show many stable conformations and 
multiple favorable ion-binding sites. 

All three binding sites, which were previously hypoth- 
esized based on the structure (Payandeh et al., 2011), 
are consistent with our simulations. In the 2-D PMF 
maps, Sin appears in the stable conformations but not 
necessarily Scen- The outer binding site S H fs is not stable 
because of the flexibility of the Glul77 side chains, but 
Glu 1 77 can form a binding site with the help of the side 
chain oxygen atom of Serl78. It can also contribute to 
the coordination shell of a Na + ion, together with the 
Leul76 side chain of the neighboring monomer in Scen- 
Finally, it can even form a binding site in partnership 
with Leu 176 of the same chain. The computational re- 
sults also indicate that the ion does not move along the 
central pore axis (Fig. 3 C) because of the relatively 
large radius of the channel (Fig. 6 A) . Instead, an ion is 
coordinated by the main chains and charged side chains 
from one or two of the monomers, and by water mole- 
cules (Figs. 3, A and C, and 6 E) . In addition, MD simula- 
tions with different ion configurations inside the channel 
pore suggest that movement of the ions is weakly coupled 
compared with the strong coupling in K + channels, which 
provide a very snug fit for their preferred ion species. 

The most important biological function of Na + chan- 
nels is to catalyze Na + permeation into the cell, while si- 
multaneously preventing permeation by K + ions. Within 
the integrated 2-D PMFs (Fig. 3 D), the barriers and wells 
depend on the occupancy in the filter. Na-Na displays 
the deepest energy wells, whereas K-K has the shallow- 
est ones in the filter, indicating that the channel selec- 
tively binds Na + . Calculation of the free energy of binding 
for two ions indicates that binding of the ion pair, Na- 
Na, is favored over K-K (Table 1). Furthermore, an en- 
tering Na + ion needs to overcome a relatively high energy 
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barrier to reach the energy minimum in the selectivity 
filter, when a K + is already present (panel Na-K) . This is 
further illustrated in the 1-D cross sections (Fig. 3 E) , 
where one ion, K + , is fixed in the filter at the position 
corresponding to the stable minimum. This example 
shows that it is energetically feasible for Na + to pass K + 
in the narrowest part of the filter, even though this 
may not represent the most probable reaction pathway 
for permeation. 

Comparison of the 2-D PMFs for E177D (Fig. 6 C) 
shows that this mutation decreased the number of po- 
tential binding sites (Fig. 6 D). Furthermore, this con- 
servative amino-acid substitution modulates the relative 
numbers of different coordinating ligands, thus altering 
local chemical moieties involved in ion coordination 
leading nearly the same relative binding free energies 
for Na + and K + to the filter. To illustrate how ion coordina- 
tion along the permeation pathway is affected by the mu- 
tation, we performed analysis of the position-dependent 
ion coordination sphere for stable binding sites found 
in 2-D PMF maps (Fig. 6 D). The average coordination 
number of Na + and K + does not change substantially 
with the E/D mutation. However, the chemical compo- 
sition of the coordination sphere does. There is a small 
decrease in the average number of carboxylates coordi- 
nating each Na + ion from the extracellular side to the 
position of Z = 2.5 A (free energy minimum for Na + ) , as 
shown in Table 2. Meanwhile, the probability of finding 
carboxylates in the potassium coordination sphere in- 
creases by ^60% in both stable sites identified by the 
PMF computations. Therefore, in the mutant, K + can gain 
additional stabilization from carboxylates, whereas this 
mode of K + coordination is largely missing in the WT. At 
the same time, K + coordination by water decreases as 
compared with the WT protein. The water ligating K + in 
the WT is replaced by carboxylates in the mutated pore. 
For example, in NaVAb E177D, the coordination with 
carboxyl oxygens was substituted by coordination with 
water molecules for Na + but not for K + (Fig. 6 D) . There- 
fore, hydration and coordination in the pore combine 
to play an important role, which appears sufficient to 
define the properties of an important ion-binding site 
and the entry barrier of bacterial Nav channels. The re- 
sults of PMF computations together with coordination 



TABLE 2 



Site-dependent ion coordination probabilities for NavAb WT and E177D 



Position of the binding site 


coo- 


—Coordination probability 


z(A) 


WT-SOD 


E/D-SOD WT-POT E/D-POT 


3.0 


0.38 


0.25 0.29 0.47 


2.5 


0.27 


0.46 0.28 0.41 



Entries represent the probabilities of finding carboxylate in the 
coordination sphere of an ion. P = 1 corresponds to all four selectivity filter 
carboxylates engaged in ion coordination. See also Fig. 6 D and its legend. 



analysis are consistent with the experimentally observed 
increase in relative potassium permeability. Thus, the 
E/D mutation allows for more favorable K + binding to 
the filter as compared with WT. There are, of course, 
limitations imposed on our analysis by the use of a closed 
or "pre-open" structure to explore the underpinnings 
of selectivity. Certainly, a prediction of absolute ion 
throughput and conductance will require detailed anal- 
ysis of an open pore structure, preferably for a channel 
from which single-channel data can be obtained. How- 
ever, by focusing our analysis on experiments evaluating 
selectivity based on ion interactions near the extracellu- 
lar entrance, we believe that we have obtained significant 
insight into the physiologically crucial discrimination 
between Na + and K + by these channels. 

Concluding comments 

The key strategy adopted by prokaryotic Nav channels 
to discriminate against K + is as follows. First, the pore is 
large enough to allow for permeation of variably hydrated 
cations. Second, within the selectivity filter, ion coordi- 
nation is highly dynamic, and a fine balance among ion 
ligation from carboxylates, carbonyls, and water mole- 
cules is at the heart of the selective entry and occupancy 
by external ions. Third, although a single binding site 
inside the filter displays a preference for Na + over K + , 
two-ion states are thermodynamically stable and are es- 
sential for rapid permeation. Finally, two procedures 
likely to modify the charge density near the narrowest 
point in the selectivity filter (E-D substitution and re- 
duction of pH) reduce the discrimination between Na + 
and K + , likely because of small changes of both energy 
minima and maxima in the PMFs for the two ions. 

Our results are consistent with a scenario in which 
differing selective pressures, during the evolution of 
prokaryotes and eukaryotes, favored a common need 
for sodium selectivity. However, each phylogenetic lin- 
eage would presumably have been subject to differing 
selection for other properties, such as details of channel 
gating kinetics, the necessity to function in different ionic 
conditions, and the presence of other complementary 
ion channels. Ultimately, the various, sometimes con- 
flicting, evolutionary demands on overlapping compo- 
nents of channel structure appear to have yielded 
different but almost equally effective "choices" of the 
molecular strategy to generate a sodium-selective con- 
ductance in prokaryotes compared with eukaryotes. 
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